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(57) Abstract: The invention relates to a computed tomography method in which an examination zone is irradiated along a helical 
trajectory by a conical radiation beam. The radiation transmitted by the examination zone is measured by means of a detector unit and 
therefrom the absorption distribution in the examination zone is reconstructed without approximations. The reconstruction comprises 
a derivation of the measuring values of parallel rays of different projections, an integration of these values along K lines, a weighting 
of these values and a baclc projection. 
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METHOD AND APPARATUS FOR EXACT CONE BEAM COMPUTED TOMOGRAPHY 

The invention relates to a computed tomography metihiod in which an 
examination zone is irradiated along a helical trajectory by a conical X-ray beam. Moreover, 
the invention relates to a computer tomograph as well as to a computer program for 
controlling the computer tomograph. 

5 In known methods of the kind set forth, utilizing approximations, the variation 

in space of the absorption or the attenuation of the radiation in the examination zone can be 
reconstructed from the measuring values acquired by a detector unit. Such approximations, 
however, give rise to artifacts in the reconstructed images, said artifacts being more 
pronounced as the angle of aperture of the radiation beam is larger in the direction of the axis 

10 of rotation ("Artifact Analysis of Approximative Cone-Beam CT Algorithms, Medical 
Physics, Vol, 29, pp. 51-64, 2002). 

Known exact methods are usually based on Radon inversion. They require a 
large amount of calculation work and give rise to discretization errors in the reconstructed 
images. 

15 Moreover, an exact method which utilizes filtered back projection is known 

from "Analysis of an Exact hiversion Algorithm for Spiral Cone-Beam CT", Physics 
Medicine and Biology, Vol. 47, pp. 2583-2597 (El). This method again requires a large 
amoimt of calculation work, thus giving rise to long reconstruction times. 

Therefore, it is an object of the present invention to provide a method which 
20 enables faster, exact reconstruction of the absorption distribution in the examination zone. 

In respect of the method this object is achieved by means of a computed 
tomography method which comprises the steps of: 

generating, using a radiation source, a conical radiation beam which traverses 
an examination zone or an object present therein, 
25 - generating a relative motion between the radiation source on the one side and 

the examination zone or the object on the other side, which relative motion comprises a 
rotation about an axis of rotation and a displacement parallel to the axis of rotation and is 
shaped as a helix. 
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acquiring measuring values which are dependent on the intensity in the 
radiation beam on the other side of the examination zone by means of a detector unit during 
the relative motions, 

reconstructing a CT image of Ifae examination zone from the measuring 
5 values, in which reconstruction an exact 3D back projection comprising the following steps is 
carried out: 

determining the partial derivative of measuring values of parallel rays with 
different radiation source positions in conformity with the angular position of the radiation 
source, 

10 - weighted integration of the derived measuring values along K lines, 

multiplying all measuring values by a weighting factor which corresponds to 
the cosine of the cone angle of the beam associated with the relevant measuring value, 

multiplying all measuring values by a weighting factor which corresponds to 
the reciprocal value of the cosine of the fan angle of the beam associated vsdth the relevant 
15 measuring value, 

reconstructing the absorption of each object point by back projection of the 
measuring values. 

In conformity witit the method which is known from El, prior to the back 
projection the measuring values must be multiplied by weighting factors which are dependent 

20 on the location of the object point to be reconstructed in the examination zone. In contrast 
therewith, in accordance with the invention the measuring values are multiplied by weighting 
factors prior to the back projection, said weighting factors being dependent on the location of 
the measuring value on the detector unit Because the number of object points to be 
reconstructed generally is much smaller than the number of detector elements, exact 

25 reconstruction is thus possible while requiring a comparatively small amoxmt of calculation 
work only. Moreover, as opposed to the method which is known from El , the integration 
interval in the back projection in accordance with the invention is not dependent on the object 
point, so that it is not necessary to determine an integration interval for each object point 
during the back projection, thus leading to a ftirther reduction of the amount of calculation 

3 0 work reqtdred. 

Claim 2 describes a preferred reconstruction method which involves an 
amoxmt of calculation work which is smaller in comparison with other methods and which 
leads to a favorable image quality. 
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Claim 3 discloses a preferred version of the integration via a K line, notably 
the multiplication by means of a weighting factor which leads to a high image quality. 

A computer tomograph for carrying out the method is disclosed in claim 4. 
Claim 5 defines a computer program for controlling a computer tomograph as claimed in 
5 claim 4. 



The invention will be described in detail hereinafter with reference to the 
drawings. Therein 

10 Fig. 1 shows a computer tomograph which is suitable for carrying out the 

method in accordance with the invention. 

Fig. 2 shows a flow chart of the method in accordance with the invention. 
Fig. 3 shows a PI line and a path of integration for a point in the examination 

zone, 

1 5 Fig. 4 shows the PI line and the path of integration for a point in the 

examination zone projected in a plane perpendicular to the axis of rotation. 
Fig. 5 shows parallel beams with different beam positions. 
Fig. 6 shows a K plane and a K line. 

Fig. 7 is a diagrammatic view of two different radiation source positions of the 
20 radiation source, a point in the examination zone and the trajectory projected in a plane 
perpendicular to the axis of rotation, and 

Fig. 8 shows the fan beams formed by rebinning in parallel planes. 



25 The computer tomograph shown in Fig. 1 comprises a gantry 1 which is 

capable of rotation about an axis of rotation 14 which extends parallel to the z direction of the 
co-ordinate system shown in Fig. 1 . To this end, the gantry 1 is driven by a motor 2 at a 
preferably constant, but adjustable angular speed. A radiation source S, for example, an X- 
ray source, is mounted on the gantry 1 . The radiation source is provided with a collimator 

30 arrangement 3 which forms a conical radiation beam 4 from the radiation generated by the 
radiation source S, that is, a radiation beam which has a finite dimension other than zero in 
the z direction as well as in a direction perpendicular thereto (that is, in a plane perpendicular 
to the axis of rotation). 
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The radiation beam 4 traverses a cylindrical examination zone 13 in which an 
object, for example, a patient on a patient table (botii not shown) or also a technical object 
may be situated. After having traversed the examination zone 13, the mdiation beam 4 is 
incident on a detector unit 16 which is attached to the gantry 1 and comprises a detector 

S surface which comprises a plurality of detector elements arranged in rows and colunms in the 
form of a matrix in the present embodiment. The detector colxmms extend parallel to the axis 
of rotation 14, The detector rows are situated in planes which extend perpendicularly to the 
axis of rotation, that is, preferably on an arc of a circle around the radiation source S. 
However, they may also have a different shape, for example, an arc of a circle aroimd the 

10 axis of rotation 14 or be rectilinear. Generally speaking, the detector rows comprise more 
detector elements (for example, 1000) than the detector columns (for example, 16). Each 
detector element struck by the radiation beam 4 produces a measuring value for a ray of the 
radiation beam 4 in each position of the radiation source. 

The angle of aperture of the radiation beam 4, denoted by the reference Omaxy 

1 S determines the diameter of the object cylinder within which the object to be examined is 
situated during the acquisition of the measuring values. The angle of aperture is defined as 
the angle enclosed by a ray situated at the edge of the radiation beam 4 in a plane 
perpendicular to the axis of rotation with respect to a plane defined by the radiation source S 
and the axis of rotation 14. The examination zone 13, or the object or the patient table, can be 

20 displaced parallel to the axis of rotation 14 or to the z axis by means of a motor 5. As an 
equivalent, however, the gantry could also be displaced in this direction. 

When a technical object is concemed instead of a patient, the object can be 
rotated during an examination while the radiation source S and the detector unit 16 remain 
stationary. 

25 When the motors 2 and S operate simultaneously, the radiation source S and 

the detector unit 16 describe a helical trajectory relative to the examination zone 13. 
However, when the motor 5 for the displacement in the direction of the axis of rotation 14 is 
stationary and the motor 2 rotates the gantry, a curcular trajectory is obtained for the radiation 
source and the detector imit 16 relative to the exanodnation zone 13. Hereinafter only tiie 

30 helical trajectory will be considered. 

The measuring values acquired by the detector unit 16 are applied to an image 
processing computer 10 which is connected to the detector unit 16, for example, via a 
wireless data transmission (not shown). The image processing computer 10 reconstructs the 
absorption distribution in the examination zone 13 and reproduces it, for example, on a 
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monitor 1 L The two motors 2 and 5, the image processing computer 10, the radiation source 
S and the transfer of the measuring values from the detector xinit 16 to the image processing 
computer 10 are controlled by a control imit 7. 

In other embodiments the acquired measuring values can first be applied to 
one or more reconstruction computers for reconstruction, said computers applying the 
reconstructed data to the image processing computer, for example, via an optical fiber cable. 

Fig. 2 shows a flow chart illustrating a version of a measuring and 
reconstruction method that can be executed by means of the computer tomograph shown in 
Fig. 1. 

After the initialization in the step 101, the gantry rotates at an angular speed 
which is constant in the present embodiment. However, it may also vary, for example, in 
dependence on the time or on the radiation source position. In the step 103 the examination 
zone, or the object or the patient table, is displaced parallel to the axis of rotation and the 
radiation of the radiation source S is switched on, so that the detector imit 16 can detect the 
radiation from a plurality of angular positions. 

In order to xmderstand the next steps, reference is made to the following 
equation from "Analysis of an Exact Inversion Algorithm for Spiral Cone-Beam CT\ 
Physics Medicine and Biology, Vol. 47, pp. 2583-2597: 

2^ |x-y(5)|.idsm/a? ^ 

This equation describes an exact reconstruction of the absorption by back 
projection of the measuring values. Therein, f(x) denotes the spatial absorption distribution in 
the examination zone in the location x and Ipj(x) describes the part of the helix which is 
enclosed by a PI line 3 1 . 

The PI line 31 of an object point 35 in the location x in the examination zone 
and /pi(,j are shown in Fig. 3 and Fig. 4 and will be described in detail hereinafter. The 

radiation source moves relative to the examination zone aroxmd an object point 35 on a 
helical path 17. The PI line 3 1 then is the line which intersects the helix in two locations and 
the object point 35, the helical segment Ipi(x) enclosed by the line then covering an angle 
smaller than 2n. 

Furthermore, in the equation (1) the reference s is the angular position of the 
radiation source S on the helix related to an arbitrary but fixed reference angular position and 
the reference y(s) is the position of the radiation source in three-dimensional space. 
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The measuring value I>f(y,®) can be described by the following line integral 



D^iyiq),®)= jdl /(y + /0). 



(2) 

The unity factor ® therein indicates the direction of the ray associated with the measuring 
value. 

5 In the step 105 in conformity with equation (1) the measuring values are 

partially derived according to q, that is, in conformity with the angular position of the 
radiation source in the location g=s. In this respect it is to he noted that only y is dependent 
on q and not ©, so that measuring values of parallel rays have to be taken into account for the 
derivation. Parallel rays have the same cone angle, the cone angle of a ray being the angle 

10 enclosed by the projection of the ray in the xz plane of the co-ordinate system shown in Fig. 
1 relative to the ray which extends through the axis of rotation and perpendicular thereto. As 
is shown in Fig. 5, rays having the same cone angle are incident on the same detector row in 
the case of a focus-centered detector, so that for the partial derivation measuring values of the 
same row but from different projections are taken into account. The derivation can then take 

1 5 place, for example, by means of the finite difference method. 

The unity factor 0 is dependent on the K angle y which can be described by 
means of the so-called K planes 5 1 . The K planes 5 1 will be described in detail hereinafter. 
In order to determine a K plane 51a function 



ms^ + (n - rn)s ^ . _ 



n 

fns + (n — m)s2 



n 



20 (3) 

is introduced, which function is dependent on non-negative, integer values n and m, where 
n>m. In this embodiment «=2 and m=L However, other values n, m may also be chosen. The 
equation (1) would nevertheless remain exact, and only the position of the K planes 51 would 
change. Furthermore, the vector function 

25 
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Ibi^i (s, s, )) - y(s)]x \yis^)- y(^)] | * 



\y{s, (s, )) - y (^)]x [y (^2 ) - vi^)] 



sgn( - J), 0 < 1^2 - j| < 271 



U(j,jj) 




(4) 



and the unity vector 



5 



p(^,x) = 



x-y(5) 
|x-y(^)l 



(5) 



are defined. The vector p then points firom the radiation source position y(s) to the position x. 
In order to determine Hie K plane, a value S2 s Ipi(x) is chosen so that y(s), y(si(s, 82)) y(s2) and 
10 X are situated in one plane. This plane is referred to as the K plane 5 1 and the line of 

intersection between the K plane 5 1 and the detector surface is referred to as the K line 53. 
Fig. 6 shows a fan-like part of a K plane. The edges of the fan meet at the location of the 
radiation source. This definition of the K plane 51 is equivalent to solution of the equation 

15 (x-y(s)).u(s,S2) = 0,S2eIpi(x) (6) 

according to S2 . Thus, u is thus the normal vector of the K plane 5 1 . In order to determine the 
vector function ©(s,x,y) the vector 



20 is defined. Using the definition for p and e, the vector function ©(s,x,y) can be expressed as 
follows: 



e(s,x) = cosy . P(s,x) + siny . e(s,x) 



(7) 



©(s,x,y) = cosy . P(s, x) + siny . e(s,x) (8) 



Because both vectors P and e are oriented perpendicularly to u, the K angle y indicates the 
direction of the vector © and hence the direction of a ray within a K plane. 



25 



The K planes and K lines are described in detail in El which is e3q)licitly 
referred to herein. 
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In the step 107 the measuring values derived along K lines are multiplied by a 
weighting factor, corresponding to the inverse sine of the K angle y, and integrated in 
conformity with the eqiiation (1). To this end, for each location x in the examination zone and 
for each projection angle there is determined a K line; as described above, a value S2 s Ipi(x) is 
5 then chosen to be such tiaat y(s), y(si(s,S2)), y(s2) and x are situated in one plane, that is, the K 
plane. The K line is then determined as the line of intersection between the K plane and the 
detector surface. The multiplications by the weightuig factor and the integrations can be 
performed, for example, by means of Fourier filtering. 

The derived and integrated measuring values can be represented by the 
1 0 following equation: 

p(y(s), (t>(s, X)) = J^^Df (y(q), ©(s,x,Y))|q.s. (9) 

Therein, p(y(s),<l>(s,x)) denote the derived and integrated measuring values and <l)(s,x) is a 
1 5 unity factor which points from the radiation source position y(s) in the direction of the 
location x in the examination zone. 

The missing integration step in the equation (1) or the back projection of the 
measuring values can now be described by the following eqxxation: 

20 In conformity with this equation each measuring value must be multiplied by 

the factor l/|x-y(s)| in order to reconstruct the spatial absorption distribution in the 
examination zone. This factor is dependent on the location x, so that it must be calculated 
anew for each combination of the radiation sovirce position y(s) and the location x. Moreover, 
the integration over s takes place along the segment of the helix Ipi(x). Thus, the integration 

25 interval is dependent on the location x in which the absorption is to be determined, so that the 
integration interval must be determined for each location x. Because the integration in 
coEdFormity with the equation (10) would require a large amoxmt of calculation work for these 
reasons, the integration variable s is replaced by the projection angle <p hereinafter. The 
projection angle 9 is then the angle enclosed by the PI line of the object point x projected in a 

30 plane perpendicular to the axis of rotation (referred to hereinafter as the xy plane) and the 
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projection on the xy plane of the ray which passes through the location x while emanating 
from the radiation source. 

The integration in the equation (10) is carried out along the helical segment 
Ipi(x). This segment is enclosed by the PI line, so that an integration must be performed for 
S each measuring value from 0 to tc after substitution of the integmtion variables. Thus, the 
integration interval is the same for each location x in the examination zone. 

The relationship between the integration variables ds and dcp can be derived 
from Fig. 7. This Figure shows a projection of a helix 17, the object point 35 in the location x 
and tihie radiation source positions y(s) and y(s4-ds) on the xy plane. The following equation 
10 results from Fig. 7: 



20 



I ^:ry(y) 1 d^cosE iJdjcose 

d<p = 



(11) 

Therein, P^y denotes the projection operator for the projection of a vector in the xy plane and 
R denotes the radius of the helix 17. The fan angle 8 is the angle enclosed by the normal from 
1 5 the radiation source position to the axis of rotation and the projection of the ray which 
emanates from the radiation source position and passes through the location x on the xy 
plane. The indices x and y describe x and y components of a vector. The components relate to 
the cartesian co-ordinate system shown in Fig. 1 . 

Thus, the following equation can be derived for the cone angle X: 



cosX=^-'''^'^^yr^\ (12) 



The equations (10), (1 1) and (12) yield 



25 m=-:^]dp-^^pb'(.smMs(q>),x)). (13) 



In conformity with this equation, in the step 109 the measuring vedues are 
multiplied by a first weighting &ctor, corresponding to the cosine of the cone angle A,, and by 
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a second weighting factor which corresponds to the reciprocal value of the cosine of the fan 
angle e. Moreover, the measuring values can be multiplied by the inverse radius R. Because 
the radius is constant during the acquisition, the latter multiplication can also be performed 
after the back projection. 

5 For small angles X and 8, the multiplication by tiie weighting factors cos(A,) 

and l/cos(8) can be ignored, because the cosine of these angles is then approximately one. 

The weighting factors in tibie step 109 are dependent on the cone angle X and 
the fan angle 8. The weighting factors are thus the same for all locations x in the examination 
zone which are traversed by the same ray, meaning that for these locations the weighting 

10 factors have to be calculated only once. In comparison with the known weighting by the 
weighting factor l/|x-y(s)| of the equation (1), the foregoing leads to a substantial reduction 
of the required amount of calculation work. 

Prior to the back projection rebinning of the measuring values can be 
performed in the step 1 1 1 . As a result of the rebinning operation the measuring values are 

1 5 resorted and re-interpolated as if they had been measured by means of a different radiation 
source (an elongate radiation source which is arranged on a part of a helix and is capable of 
emitting each time mutually parallel fan beams) and by means of a different detector (a flat, 
rectangular "virtual" detector containing the axis of rotation 14), 

This will be described in detail with reference to Fig. 8. Therein, tihie reference 

20 numeral 17 denotes the helical trajectory wherefrom the radiation soxirce irradiates the 
examination zone. The reference numeral 43 denotes a fan-shaped radiation beam which 
emanates from the radiation source position So and whose rays propagate in a plane 
containing the axis of rotation 14. The conical radiation beam emitted by the radiation source 
in the position So can be assumed to consist of a plurality of flat fan beams which are situated 

25 in planes which are parallel to the axis of rotation 14 and intersect in the radiation source 
position So. Fig. 8 shows only a single one of these fan beams, that is, the fan beam 43. 

Moreover, Fig. 6 shows further fan beams 41, 42 and 44, 45 which extend 
parallel to the fan beam 43 and are situated in planes which are parallel to one another and to 
the axis of rotation 14. The associated radiation source positions S.2, S.i and Si, S2 are 

30 occupied by the radiation source S before and after having reached the radiation source 

position So, respectively. All rays in the fan beams 41 to 45 have the same projection angle. 

The fan beams 41 to 45 define a radiation beam 70 having a tent-like shape. 
For each group of fan beams there is defined a rectangular, virtual detector 160 which is 
situated in a plane which contains the axis of rotation 14 and is oriented perpendicularly to 
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the parallel fan beams of a group. The comer points of the virtual detector 160 constitute the 
pxmcture points of the rays, incident on the oppositely situated helical segment from the outer 
radiation source positions, through this plane. For the fan beam 70 in Fig. 8 these are the 
points of intersection of the fan beams 41 and 45 with the helix. On the rectangular detector 
5 160 detector there are defined elements whicK are arranged in a cartesian fashion, that is, in 
rows and columns, on which the measuring values are re-inteipolated. 

The measuring values obtained after the rebinning are subsequently used for 
the reconstruction of the absotption distribution in the exanoination zone by a back projection 
which is in this embodiment in conformity with the equation (13). 

10 In the step 113 a voxel V(x,y,z) is determined within a selectable (x,y,z) zone 

(field of view or FOV). Subsequently, in the step 115 a projection angle is selected within the 
range [cpo, cpo + ti] is selected, where 90 is the angle at which the voxel V(x,y,z) enters the 
radiation beam. In the step 1 17 it is checked whether a ray of the projection extends through 
the center of the voxel V(x,y,z). If no ray of the projection passes through the center of the 

1 5 voxel, the associated value must be deteraxined by interpolation of the measuring values of 
neighboring rays. The measuring value that can be associated with the ray passing through 
the voxel, or the measuring value obtained by interpolation, is accumulated on the voxel 
V(x,y,z) in the step 1 19. In the step 121 it is checked whether all projections with the 
projection angles 90 to <pO-Hc have been taken into account. If this is not the case, the flow 

20 chart branches to the step 115. Otherwise, it is checked in the step 123 whether all voxels 

V(x,y,z) in the FOV have been dealt with. If this is not the case, the procedure continues with 
the step 113. However, when all voxels V(x,y,z) in tihe FOV have been dealt with, the 
absorption has been determined in the entire FOV and the reconstruction method is 
terminated (step 125). 
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CLAIMS: 



1 . A computed tomography method which comprises the steps of: 
generating, using a radiation source (S), a conical radiation beam (4) which 

traverses an examination zone (13) or an object present therein, 

generating a relative motion between the radiation source (S) on the one side 
5 and the examination zone (13) or the object on the other side, which relative motion 

comprises a rotation about an axis of rotation (14) and a displacement parallel to the axis of 
rotation (14) and is shaped as a helix (17), 

acquiring measuring values which are dependent on the intensity in the 
radiation beam on the other side of the examination zone by means of a detector unit (16) 
1 0 during the relative motions, 

reconstructing a CT image of the examination zone (13) from the measuring 
values, in which reconstruction an exact 3D back projection comprising the following steps is 
carried out: 

determining the partial derivative of measuring values of parallel rays with 
1 5 different radiation source positions in conformity with the angidar position of the radiation 
source, 

weighted integration of the derived measuring values along K lines, 
multiplying all measuring values by a weighting factor which corresponds to 

the cosine of the cone angle of the ray associated with the relevant measuring value, 
20 - multiplying £dl measuring values by a weighting factor which corresponds to 

the reciprocal value of the cosine of the fan angle of the beam associated with the relevant 

measuring value, 

reconstructing the absorption of each object point by back projection of the 
measuring values. 

25 

2. A computed tomography method as claimed in claim 1, in which in the 
reconstruction step rebinning of the measuring values is performed prior to the back 
projection so as to form a number of groups, each group comprising a plvirality of planes 
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which extend parallel to one another and to the axis of rotation and in which a respective fan 
beam (41 ... 45) is situated. 

3 . A computed tomography method as claimed in claim 1 , in which the 
3 integration of the measuring values along K lines comprises the following steps: 

determining a K plane for each radiation source position and each location to 
be reconstructed in the examination zone, 

determining K lines, that is, lines of intersection between the K planes and a 
detector surface of the detector unit, 
10 - multiplying the measuring values on each K line by a weighting factor which 

corresponds to the reciprocal value of the sine of a K angle, 

integrating the measuring values along the K lines. 

4. A computer tomograph for carrying out Ihe method claimed in claim 1, 
15 comprising: 

a radiation sotirce (S) and a diaphragm arrangement (3) which is situated 
between the examination zone (13) and the radiation source (S) in order to generate a 
radiation beam (4) which traverses an examination zone (13) or an object present tiierein, 
a detector xmit (16) which is coupled to the radiation source (S), 
20 - a drive arrangement (2, 5) which serves to displace an object present in the 

examination zone (13) and the radiation sotirce (S) relative to one another about an axis of 
rotation (14) and/or parallel to the axis of rotation (14), 

a reconstruction xmit (10) for reconstructing the spatial distribution of the 
absorption within the examination zone from the measuring values acquired by the detector 
25 unit (16), 

a control unit (7) for controlling the radiation source (S), the detector unit (16), 
the drive arrangement (2, 5) and the reconstruction unit (10) in conformity with the steps 
disclosed in claim 1 . 

30 5. A computer program for a control unit (7) for controlling a radiation source 

(S), a diaphragm arrangement (3), a detector unit (16), a drive arrangement (2, 5) and a 
reconstruction unit (10) of a computer tomograph so as to execute the steps disclosed in 
claim 1. 
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